Non-afRne deformation of inherent structure as static signature of cooperativity in 

supercooled liquids 
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We unveil the existence of non-afhnely rearranging regions in the inherent structures (IS) of 
supercooled liquids by numerical simulations of two- and three-dimensional model glass formers 
subject to static shear deformations combined with local energy minimizations. In the liquid state 
IS, we find a broad distribution of rather large rearrangements which are correlated only over 
small distances. At low temperatures, the onset of the cooperative dynamics corresponds to much 
smaller displacements correlated over larger distances. This finding indicates the presence of non- 
affinely rearranging domains of relevant size in the IS deformation, which can be seen as the static 
counterpart of the cooperatively rearranging regions in the dynamics. This idea provides new insight 
into possible structural signatures of slow cooperative dynamics of supercooled liquids and supports 
the connections with elastic heterogeneities found in amorphous solids. 

PACS numbers; 61.43.Fs,64.70.Q-,05.20.Jj 



When a liquid is cooled down to its glass transition 
temperature, particle motion slows down enormously and 
becomes highly cooperative [ij: The more the system 
becomes glassy, the more the relaxation requires coop- 
erative rearrangements of a large number of particles. 
Experimental and theoretical investigations as well as 
simulation results on dynamical heterogeneities strongly 
support the presence of cooperatively rearranging regions 
of growing size P, H, d, Q • However, it is still debated 
whether and how the onset of such cooperative dynamics 
can take place without any apparent or straightforward 
connection to structural changes. Such links have been 
searched for in various directions. For example, different 
connections between the local propensity of particles to 
motions and the underlying structural features of the sys- 
tem have been investigated, but they still remain quite 
elusive 

An alternative (and so far complementary) description 
of the dynamics of supercooled liquids is based on the 
way the system explores its potential energy surface @. 
When the system enters the supercooled regime, the pres- 
ence of local minima, or inherent structures (IS), turns 
out to be extremely important and to produce complex 
trajectories in the potential energy landscape: The time 
evolution of the system samples the basins associated to 
the local minima, or groups of such basins, on a rela- 
tively short time scale, whereas transitions between dif- 
ferent groups of basins separated by major energy bar- 
riers take a much longer time and imply a slow, aging 
process, during which the system is out-of-equilibrium. 
Such a description, however, does not fully account for 
the cooperative dynamics associated to dynamical het- 
erogeneities. Yet, possible connections have been dis- 
cussed Jyland recent works have gathered additional in- 
sights [fi,i,0,[ni,[ii. 

In the present paper, we report about a signature of the 
onset of cooperative dynamics in the structural features 
of supercooled liquids from a novel perspective. The ap- 
proach is motivated by a recent theory [l3| - based on a 



general framework of nonequilibrium thermodynamics - 
which relates the response of a glassy system to an ap- 
plied deformation to the corresponding change of its IS. 
In contrast to previous approaches, here it is suggested 
that the reversible part of the dynamics changes con- 
siderably when approaching the glass transition. Since 
the reversible dynamics basically mirrors the space trans- 
formation behavior of the system, the suggestion in 
guides the way toward investigations of static rather than 
dynamical properties. Accordingly, the most important 
feature of IS transformation is a hindrance to affine de- 
formation: The typical linear size of the regions where 
non-affinc deformations take place is roughly the particle 
size in the liquid and is expected to become larger upon 
approaching the glass transition. 

We explore these theoretical considerations by com- 
puter simulations of different model systems for super- 
cooled liquids subject to static shear deformations with 
small amplitudes. In the liquid phase, a broad distri- 
bution of rearrangements is found, centered around a 
rather large mean value. When the system enters the 
supercooled regime, instead, its IS displays the onset of 
an enhanced collective behavior with much smaller dis- 
placements. This feature, which does not depend on the 
particular model or deformation considered, has been de- 
tected here for the first time and indicates the existence 
of non-afhnely rearranging regions of relevant size in the 
IS of supercooled liquids. The presence of such collective 
rearrangements in the IS response is strongly evocative 
of the cooperativity characterizing the dynamics, and it 
is in fact observed in the same range of temperatures. On 
this basis, we propose that these non-afhnely rearranging 
regions are the IS counterpart of the cooperatively rear- 
ranging regions observed so far only in the dynamics. 

Methods and numerical simulations. We consider three 
well established model systems for supercooled liquids: 
(a) The three-dimensional (3D) binary Lennard- Jones 
mixture of Ref . [l^ ; (&) The 3D soft spheres binary mix- 
ture studied in Ref. (c) The 2D soft spheres binary 
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mixture in Ref. [16[. In all three models, the densities 
have been chosen such that crystallization is prevented. 
Upon cooling, these systems therefore go from the liquid 
to the supercooled and glassy regime. In order to inves- 
tigate the response of IS to an applied deformation we 
have designed the following procedure First, start- 
ing from a particle configuration X = {vi \ equilibrated at 
a given temperature T, wc deform it afhncly — > and 
subsequently find the inherent structure of the deformed 
configuration X^'^ = {r^''}. For the energy minimization 
we use a conjugate gradient algorithm. Second, starting 
from the same initial configuration X, we first find the 
IS corresponding to the initial configuration X^ = {r^} 
and subsequently apply the same afRne deformation to 
obtain X^'^ = {r^'*}. All subsequent analysis is based on 
the comparison between the two configurations X'^'^ and 
X^"^, which is quantified by the mismatch, i.e. the non- 
affine displacement field {di = r^'^ — r^'^}. In the limit 
of a flat energy landscape, {dj} is identically zero. By 
lowering the temperature into the supercooled regime, 
we expect the presence of basins and barriers in the en- 
ergy landscape to introduce collective contributions to 
the displacement field. We carefully chose the deforma- 
tion amplitudes sufficiently large to make the system typ- 
ically leave a basin. Supplementary information is given 
in the EPAPS document We have considered two 

different kinds of deformation: (i) The spatially modu- 
lated shear deformation ^ = -f 7sin(A:„yi)ej,, 
where 7 denotes the maximum amplitude of the defor- 
mation and kn = wk / L the wave- vector with L the 
size of the simulation box; [ii) The homogeneous shear 
Ti ^ rf = Yi -\~ "fUiex, where Lees-Edwards boundary 
conditions are used. We have applied the deformations 
(i) and (ii) to the systems (a), (6) and (c) with 7 varying 
between 10~^ and 10°. All quantities are given in re- 
duced Lennard-Jones units. For all the systems we have 
prepared 20 to 50 independent samples of N particles 
(N varying from 500 to 8000) which have been care- 
fully equilibrated AH quantities discussed are aver- 
aged over the independently prepared samples and the 
error bars obtained from sample-to-sample fluctuations. 
Some simulations were performed using the LAMMPS 
code [l^. The following analysis concerns temperatures 
T in the range from the high T Arrhenius dependence of 
dynamic quantities to the onset of caging and cooperative 
dynamics. 

Results. A first, qualitative understanding can be ob- 
tained via a simple visualization of the mismatch field 
{dj} at different temperatures and for different ampli- 
tudes of deformation. In Fig. [H {di} vectors arc plotted 
for the system (c) and deformation (ii) at high tempera- 
ture T = 1.0 (left) and low temperature T = 0.46 (right). 
The length of the vectors are suitably enlarged for clar- 
ity. At high T, for a fixed amplitude of deformation, the 
mismatch vectors display a relatively large magnitude as 
well as uncorrelated directions. This indicates that, at 
these temperatures, the local changes to the IS induced 
by our procedure are in fact dominated by the thermal 




FIG. 1: Mismatch field {di} for the 2D soft spheres binary 
mixture of Ref. ^ at T = 1.0 (left) and T = 0.46 (right) for 
homogeneous shear deformation with amplitude 7 — 0.001. 
For better visibility, the lengths of the vectors are multiplied 
by a factor 14 and 600 for T = 1.0 and T = 0.46, respectively. 



fluctuations of the system. At low temperatures, instead, 
{di} have much smaller magnitudes for the same 7 (see 
figure caption for details) but, at the same time, they 
appear rather spatially correlated in magnitude and in 
direction. These two kinds of behavior are ubiquitous in 
our numerical study and appear systematically for a wide 
range of 7 values [l8|. Interestingly enough, in all cases 
considered, the transition from the one to the other be- 
havior is coupled to the onset of the supercooled regime 
(for example, for the system (c) in Ref. [TB| the onset of 
the glassy dynamics can be set around T w 0.46). These 
findings support the idea that a qualitative change of the 
response of IS to deformation is strongly correlated to 
the onset of the cooperative dynamics [7|. 

A first, simple characterization is given in terms 
of the mean length of the mismatch vectors Id = 
E^ df In Fig.H Uj, T)//rf(7, T = 1), obtained 
in the system (&) using the deformation (i) of amplitude 
7 and wave vector ki = tt/L, is plotted as a function of 
the temperature T. For sufficiently small values of 7, the 
relative global magnitude of the mismatch field shows a 
rather steep decrease with decreasing temperature. In 
fact in Ref. , it has also been observed that the critical 
deformation amplitude to cause IS transitions strongly 
depends on temperature. Remarkably, there are clearly 
two well distinguished temperature regimes for the typi- 
cal length of the mismatch vectors. In addition, the tran- 
sition from the high temperature to the low temperature 
one becomes steeper for small 7 and takes place at the on- 
set of the cooperative dynamic regime {T ~ 0.4), i.e. very 
close to Teis, where the IS energy starts to strongly de- 
pend on temperature, and well above structural arrest 
(Tmct^ 0.303) [ll,[l3. 

An estimate of the fraction of particles involved in the 
rearrangements described by the mismatch field can be 
obtained by the dimensionless width of the length dis- 
tribution n = (7V-i[E,df]2/E4df]2) (the participation 
ratio of Ref. [l^). In the main frame of Fig. O If is plot- 
ted as a function of temperature T for 7 = 10~^, 10"'^ 
and 10~^ for the same system and the same wave vector 
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FIG. 2: (color online). Normalized mean length of mismatch 
field ld{y,T)/ldij,T — 1) as a function of T in the 3D binary 
mixture of Ref. [1^, obtained with the spatially modulated 
deformation ( i) of amplitude 7 and wave vector ki . 



as in Fig. [21 The width 11 increases with decreasing tem- 
perature below T — 0.4. In particular, for the very small 
values of 7 shown in the figure, 11 displays a rather sharp 
increase and shows that the fraction of particles involved 
in such non-affine deformation of the IS increases from 
10% to 40% — 60%. In order to further characterize the 
nature of the non-afhne rearrangements we have mea- 
sured the average degree of correlation C(r) between the 
directions of the mismatch vectors of different particles i 
and j separated by a distance r, C(r) = (d^ •dj)/(d^). In 
the inset of Fig.[3l C(r), obtained for 7 = 0.01 in the sys- 
tem (6) with the deformation (z), is plotted as a function 
of r. The data show that the low temperature regime 
is characterized by a definitely larger degree of correla- 
tion at the same distance r. Moreover, the distribution 
of values of d^ ■ dj/{df) for pairs of particles separated 
by distance r ± dr, plotted in Fig. [H widens strongly 
at low temperature, indicating that the extended corre- 
lation in the response of IS to deformation corresponds 
to a high degree of heterogeneity. This feature is strik- 
ingly similar to what is observed in the dynamics with 
the average time correlation and the fluctuations around 
this average, generally used to quantify dynamical het- 
erogeneities. But it is observed here, in the same range 
of temperatures (l7j . in a static quantity. On the whole, 
the emerging new scenario is that the onset of coopcra- 
tivity in the real-time dynamics corresponds to growth of 
correlated and heterogeneous domains in the non-affine 
deformations of the IS. 

The same qualitative behavior is obtained in the entire 
range 10^^ < 7 ^ 0.1 and in the systems (a) and (c) 
for both deformations (i) and {ii). We have also found 
in all the systems considered that, when 7 > 1.0, the 
mismatch field always closely resembles the one obtained 
at high temperature. This indicates that beyond 7 2± 1.0 
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FIG. 3; (color online) . Mainframe; The participation ratio 11 
as a function of the temperature for the 3D binary mixture of 
Ref. p^ . in the case of the spatially modulated deformation 
(i) of amplitude 7 = 10"'*, 10~^, 10~^. Inset: Average degree 
of correlation C(r) of the direction of the mismatch vectors 
separated by a distance r for 7 = 0.01. 
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FIG. 4: (color online). Distributions of degree of correlation 
C(r) between particles at a distance r = 5 in particle diame- 
ter units. The same system and the same spatially modulated 
deformation (i) as in Fig. [3] was chosen with deformation am- 
plitude 7 = 0.01. 



the amplitude of the deformation is sufficient to erase 
any distinctive feature of the low temperature structure 
of the IS, as in a sort of Lindemann criterion Finally, 
the details of temperature dependence of the different 
quantities, e.g. the crossover width or the values of 11, 
depend to some extent on the specific system considered 
(e.g. the crossover is always smoother in the 2D system) 
and on the type of deformation [l^ . 
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Conclusions. The presented results reveal for the first 
time that the onset of the glassy dynamics corresponds to 
the onset of cooperative heterogeneous rearrangements 
in the response of IS to external, static deformations. 
In fact we have shown that in the liquid state, as a re- 
sponse to the external deformation, the IS displays large 
and uncorrelated rearrangements, dominated by thermal 
fluctuations. Near the onset of the cooperative dynam- 
ics, the response of the IS to small external deformation is 
instead characterized by small, collective non-affine rear- 
rangements involving a relevant fraction of particles . 
To our knowledge, this feature of the IS had not been rec- 
ognized until now. It has been possible here, thanks to 
the procedure designed in Rcf. [iJl which strongly relies 
on the many-body nature of the IS construction. An in- 
teresting advantage of our approach is in avoiding calcu- 
lations based on the dynamics, which usually require very 
long simulation times. The extent of such an advantage 
can be better quantified by systematically comparing our 
procedure to the dynamics [18[ . 

The heterogeneous domains of the IS, determined by 
collective non-affine displacements, are strongly evoca- 
tive of the cooperatively rearranging regions observed in 
the dynamics and are in fact detected in the same range 
of temperatures. We propose therefore that they are ac- 
tually the IS counterpart of the cooperatively rearrang- 
ing regions. If this is the case, they would deliver a static 
correlation length which grows significantly in the su- 



percooled regime, where the dynamic correlation length 
related to dynamic cooperativity starts to grow. From 
this point of view, this study suggests a new path for 
the investigations of structural signatures of glassy dy- 
namics. A more sophisticated analysis of the non-affine 
displacement field, for example, adds useful insights in 
this direction . Finally, our results have also interest- 
ing connections with the studies of elastic heterogeneities 
in amorphous solids 0, M, HI, HI, H [H, Hi, [13 . Their 
typical length scale has actually been related to the non- 
affine part of the displacement field gen erated in response 
to an elastic deformation. In Ref. [26|, where a quanti- 
tative estimate of this length scale is obtained in a real- 
istic model of amorphous silica, the connection between 
the typical length scale of elastic heterogeneities in the 
amorphous solid and the length scale typical of dynamical 
heterogeneities (which is actually of comparable magni- 
tude) in the supercooled liquid is suggested. Our results 
support this idea and offer a way to extend those studies 
to finite temperatures. On approaching the supercooled 
regime from the high temperature fiuid, we find a sig- 
nature of the onset of dynamic cooperativity in a purely 
static quantity: the non-affine, collective rearrangements 
of the IS show a striking similarity to cooperative, dy- 
namical phenomena. 
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APPENDIX A: DETAILS OF NUMERICAL 
SIMULATIONS 

For the two 3D systems (a) and (&). we have performed 
NVE and NVT molecular dynamics simulations of sys- 
tems containing = 2048 particles. For system (a) and 
(&), the chosen number densities and correspon ding size 
of the simulation b ox are p = 1.2, L = 11.95 [Alj and 
p = 0.6, L = 15.06 |A2l |. respectively. For the 2D system 
(c), wc have performed constrained molecular dynamics 
simulations in the iVPT-ensemble (with N = 500, 2048, 
and 8000), with Nose-Hoover thermo- and barostats en- 
suring on average a fixed t emp erature (0.3 < T < 2) and 
pressure {P = 13.5) value [A3j |. 
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APPENDIX B: THE MISMATCH DISTRIBUTION 

The motivation for studying the mismatch field is the 
idea, presente d wi thin the thermodynamic modeling ap- 
proach of Ref. [A5| | , that the behavior of the system under 
deformation can be described using the underlying inher- 
ent structure. Since we here study the liquid regime, we 
have chosen the mismatch field {d ~ X'^'^ — X'^'^} in order 
to subtract the affine part of the deformation, such that 
d only accounts for non-affinc contributions to the de- 
formation. For a crystalline solid, instead, the mismatch 
has to be more conveniently considered with respect to 
X^, i.e. the IS of the original configuration. In this case, 
d is therefore a measure of the difference between two 
configurations of the solid. 

There are number of related and special cases which 
help to clarify the meaning of the mismatch vectors {dj}, 
or equally, the mismatch vector field d(r): (A) For the 
special case of a perfectly flat energy landscape, d van- 
ishes identically; (B) In a real fluid, the mismatch vector 
field d does not vanish even at high temperatures, in fact 
its magnitude tends to become large in that case, but 
it will be orientationally correlated only over short dis- 
tances and vanish when averaged over volume elements 
of a few particle diameters; (C) In the case of extremely 
small deformations that do not allow to overcome even 
the slightest energy barriers, d(r) is identical to minus 
the imposed deformation at location r. This limit has 
not been reached for the range of deformation amplitudes 
7 that we have considered (see also next section); (D) 
Applying rather large deformations 7 > 10~^ within our 
procedure (but out of focus of the present study) gives in- 
formation abo ut si zes of inherent structure basins, similar 
to the work of [A6| : (E) For very low temperatures, where 
the initial configuration is an amorphous solid, X X'^, 
the mis mat ch field coincides with the quantity studied 
in Ref. [A7l |. where it was used in order to extract the 
non-affine displacements in amorphous solids. 

From the measured mismatch fields {dj} we have fur- 
ther compiled ensemble averaged probability distribu- 
tions of mismatch strength; denoted as h{d), where d 
is the norm of the mismatch vector. For case (C) con- 



FIG. 5: Scaled probability distribution function of mismatch 
strengths, with q = 2/{'yL), for T = 0.46 and 7 = 10"^ for 
system (c). Also shown is the reference distribution one would 
observe for extremely small deformations, as described in the 
text part. 
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FIG. 6: Same as Fig. [5]for T = 0.8 and 7 = 10" 



sidered above, the distribution is flat and vanishes for 
d > "fL/2 (for the case of homogeneous shear strain ap- 
plied of size 7 to a homogeneous system of linear size 
L). In Figs. [5] and m we quantitatively demonstrate for 
two representative temperatures and strain amplitudes 7, 
that even for deformations as small as those required for 
the approach presented in the manuscript, this regime 
is not observed. Moreover, the figures show that the 
distributions are characterized by an 'initial' (at small 
7) parabolic increase with d, a 'later' decrease well de- 
scribed by a power law, and a finite range of displacement 
strengths which ensures intcgrability and calculation of 
the moments 11 and presented in Figs. 2 and 3 of the 
manuscript. The occurrence of a significant fraction of 
values of qd between 1 and 10 in Fig. 6 underlines the 
enormous freedom of rearranging individual particles in 
the liquid state. 

The strong discrepancy between measured distribu- 
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tions and the reference distribution, as shown in Figs. [5] 
and [B] for two different temperatures, sheds some Hght 
on the fact that a picture of energy valleys being large 
enough to accommodate small deformations without af- 
fecting the inherent structure begins to fail at deforma- 
tions small compared to the ones needed to recognize the 
occurrence of the glass transition. 



indications, we have calculated conditional probability 
distributions for the propensity and mismatch lengths. 
We have also evaluated, at different temperatures, differ- 
ent correlation coefficients between propensity and mis- 
match length distributions, the latter calculated for dif- 
ferent values of deformation amplitudes. Wc can rule 



APPENDIX C: PROPENSITY OF PARTICLES 
TO MOTION 

For the model of Ref. jA2j | , to which most of the data 
shown in the paper refer, no study of dynamical het- 
erogeneity was available in the literature. Wc have cal- 
culated the local propensity to motion p as defined in 
Ref. jA9l | from the distribution of particle displacements 
after a time interval ~ l.Srg, where Te is the time at 
which the intermediate scattering funcion, measured at 
the wave vector of the first peak in the structure factor, 
has decayed to 1/e. The distribution is calculated from 
Nr = 20 separate simulation runs over a fixed time in- 
terval, all starting from the same particle configuration 
but with momenta randomly assigned from the Maxwell- 
Boltz man n distribution at a certain temperature. In 
Ref. |A9l | it is shown that the probability distribution 
f{p) becomes asymmetric and much wider at the onset 
of the supercooled regime. In Fig. 7 we have plotted the 
distribution of propensity values at various temperatures: 
the data show the change in the shape of the distribution 
at temperatures below T — 0.4 and higher than Tmct- 
The distribution of degree of spatial correlation plotted 
in Fig. 4 of the manuscript displays a strikingly similar 
qualitative change in the same range of temperature. 

APPENDIX D: SEARCH FOR CORRELATIONS 
BETWEEN DISTINCT APPROACHES 

One of the possible developments of this work is in the 
direction of quantitatively characterizing correlations be- 
tween the cooperative domains observed in the response 
of IS and the heterogeneities observed in particle dynam- 
ics. In principle, this would require a systematic analy- 
sis of the different measur es of dynamical heterogeneities 
proposed in the literature |A8l ]. but, to get started, let us 
consider the local propensity to motion p discussed in the 
previous secion. A significant change of shape is also ob- 
served for the distribution h{d) of magnitudes of the mis- 
match vectors, as discussed in the previous section, and 
even more strikingly for the distribution of (d^ -dj), /(C), 
see Fig. 4 of the manuscript. Starting from these first 




FIG. 7: Distribution of particle displacements (propensity to 
motion) at time r ~ 1.5re at different temperatures. 



out a clear signal of a simple and straightforward cor- 
relation between \dj\ and pj. This is illustrated e.g. in 
Fig. [SI where the propensities and mismatch vectors are 
visualized on top of the particle configuration for the 2D 
system ( c) . Correlations between particle propensity and 
the spatial correlation of the mismatch vectors give in- 
stead more encouraging indications, as also suggested by 
the cooperative nature of the mismatch vectors at low 
temperature. By defining the coarse-grained mismatch 
field D{r) = (d^(r))j, where indicates the average 
over the Nj p articles contained in the volume Vj of lin- 
ear size r [A7l | , and the analogue coarse-grained propen- 
sity field P{r), we have calculated their spatial correla- 
tion {D{r)P{r)), plotted in Fig. 8 for the model system 
(b). Such correlations apparently become longer ranged 
at low temperatures. The same behavior is observed in 
different model systems, and for different shear ampli- 
tudes. 
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FIG. 8: Visualization of particle propensity and mismatch vectors for the system (c). Increasing propensity is color coded from 
violet (dark) to red (bright). Same conditions are chosen as in Fig. 1 of the manuscript with T = 0.4 (left) and T = 1.0 (right). 
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FIG. 9: Spatial correlations between the coarse grained fields 
D{r) and P{r) as a function of the coarse graining distance 



